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ABSTRACT 



We argue that the asymmetric morphology of the blue and red shifted components of the 
outflow at hundreds of AU from the massive binary system rj Carinae can be understood from 
the collision of the primary stellar wind with the slowly expanding dense equatorial gas. Recent 
high spatial observations of some forbidden lines, e.g. [Fe III] A4659, reveal the outflowing gas 
within about one arcsecond (2300 AU) from rj Car. The distribution of the blue and red shifted 
components are not symmetric about the center, and they are quite different from each other. 
The morphologies of the blue and red shifted components correlate with the location of dense 
slowly moving equatorial gas (termed the Weigelt blob environment; WBE), that is thought to 
have been ejected during the 1887 - 1895 Lesser Eruption (LE). In our model the division to the 
blue and red shifted components is caused by the postshock flow of the primary wind on the two 
sides of the equatorial plane after it collides with the WBE. The fast wind from the secondary 
star plays no role in our model for these components, and it is the freely expanding primary 
wind that collides with the WBE. Because the line of sight is inclined to the binary axis, the 
two components are not symmetric. We show that the postshock gas can also account for the 
observed intensity in the [Fe III] A4659 line. 

Subject headings: stars: mass loss — stars: winds, outflows — stars: variables: other — stars: 
individual (Eta Car) 



INTRODUCTION 



The colliding- wind binary system rj Car is a special binary system in our Galaxy (e.g., Davidson fc Humphrevj 

19971 ). It is composed of a primary Luminous Blue Variable (LBV), a very massive star (Mi > 120Mq) 
whose on the edge of instability, and a somewhat evolved massive 0-star secondary (M2 > 3OM0). The 
binary system has a complicated temporal behavior and a complex circumbinary nebula. Both are con- 
nected to the massive and energetic winds of the system at present and in the past. Interesting dramatic 
periodic changes, termed collectively as "the spectroscopic event", take place every ^ 5.54 years close to 
periastron passage. The effect of the secondary gravity when the secondary approaches the primary on its 
highly eccentric (e > 0.9) orbit is thought to be strongly related to the occasion of the spectroscopic event. 
During the s pectroscopic event ma ny emission and absorption lines and bands show considerable changes for 
a few weeks ( Damineli et al. 20081 and references therein), e.g., a deep minimum in th e X-ray emission. The 



X-ray emission, for example, comes from the colliding stellar winds, ( Corcoran 
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2008t IPittard fc Corcoranl [iooj llQQSt iBehar et all boOTt iTeodoro et all |2012| ). while its minimum is at- 
tributed to the suppression of the secondary wind near periastron passages. This suppression occurs when 
the colliding-wind region comes very close to the secondary star a s the two stars approa ch each other, and 
the secondary gravity acts on the dense postshock primary wind teashi fc Sokeil lioOQbl). This accretion is 
aided by the formation of clumps in the primary wind (jAkashi et al.ll2011l ). Another explanation (which 
requires our less favor able orientation; see below) is that the colliding-wind region is occulted by the primary 
( Corcoran et al. 201Cll). Howev er, occulation alone cannot account for the X-ray suppression ( Kashi fc Soker 



2n09a 



Hamaeuchi et al.ll2007h . 



Apart from the colliding-wind region there is a more extended structure around the binary system. 
The complexity of the structure poses a challenge in understanding observations from the system, and 
how different parts of it were formed. The bipolar Homunculus nebula surrounding the system is the most 
prominent structural feature. It was formed during the 183 7.9 ^ 1858 Great Eruption (GE) when the primary 
LBV ejected 12 - iOMp,, a considerable part of its mass JComez et ahl bood boid ISmith fc Owockilbood 
Smith et al.ll2003llSmithfc FerlandlboOTt iKashi fc Sokeilboiol ). Both the GE and the weaker 1887.3 - 1895.3 



Lesser Eruption (LE) that fol lowed it are the source of most of the circumbinary medium. The LE was a 
much less en ergetic eruption ( Humphrevs et al. 19991 ) and only 0.1 - IMq were ejected from the primary 
(|Smithlb005l ). Part of this material was ejected inhomogeneously in equatorial directions. Material from the 
GE in the equatorial direction created the equatorial skirt, while material ejected from t he LE is thought t o 



be the source of dense blobs and the gas in their sur roundings closer to the binary svstem ([Smith et al 



These blobs are known as the Weigelt blobs (WBs; IWeigelt fc Ebersbergeiill986l : iHofmann fc Weigelt 



2004 ) 



1988) 



The WBs are concentrated on the NW side of p Car, but there is also dense material in other direction s 
fe.g.- ISmith fc Gehrjiiooot lOorland et allbooi Ichesneau et"al]b005l: ICuU et aLlliooi iMehner et allboiol) . 



In this paper we focus on the interaction of the primary wind with material in the Weigelt blob en- 
vironment (hereafter WBE). Observations of ionized iron lines from the system in recent years provide us 
with detailed inform ation about the WB E. From HST/STIS spectra of Fe II, [Fe II], [Ni II] and [Ni III] lines 
from WBs C and D . ISmith et al.l (|2004r) found the radial velocity of the WBs to be 40 km s~^. From that 



Smith et al.l (|2004l ) concluded that the WBs were originated in the LE. On the other hand, iDorland et al 



( 20041 ) presented HST observations of WBs C and D together with simulations of their propagation and 
concluded that they formed during the period 1910 - 1942, preferably during a brightening event which took 
place in 1941. For the purpose of this paper it is immaterial when the dense WBE was ejected. 

GuU et aD toO^ presented observations of broad (^ 500 km s ^) [Fe II] emission line structures which 



extend to ~ 1600 AU from the binary system, and [Fe III], [Ar III], [Ne III] and [S III] lines which extend 
up to a projected distance of only ^ 700 AU from NE to SW. The latter lines showed radial velocities 

-500 to +200 km s"^^ 



with tendency to the blue (^ 



Gull et al.l (|2009f ) observed all those forbidden lines 



disappearing during the spectroscopic event. Madura et al. |2012l ) took HST/STIS spectra of [Fe HI 



emission 

lines from slits in different position angles through the central source of ry Car. The observations covered the 
2003.5 sp ectroscopic event of r ) Car. While the low velocity component of these lines was associated with 
the WBs, iMadura et al.l (|2012r) interpreted th e high velocity co mponents of these lines as being formed in 
the wind collision zone of the binary system. iGull et al.l (|201l[ ) presented more HST/STIS high- ionization 
forbidden-line observations taken after the 2009 event, this time in the form of complete intensity maps. 
The high velocity components were evident as well. Based on comparison with SPH simulations and under 
the assumption ( w hich we will d i spute b elow) that the hues origin ate in the b inary coUiding-wind region. 



Gull et all (|2009l) . iMadura et al.l (|2011l ). iMadura et al.l (|2012D and iGull et al.l (|201ll) . concluded that the 



high velocity components of the forbidden-line observations fit an orientation in which the secondary star is 
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in the direction of the observer for most of the binary orbit. Namely, the longitude angle (the argument of 
periapsis) is w ~ +270° according to their interpretation. 



There are several other papers claiming a longitude angle of w ~ +240° - +270° (e.g., Moffat & Corcoran 


2009; Okazaki et al. 2008; Henlev et al. 2008; Hamaeuchi et al. 


2007; Nielsen et al. 2007; Gull et al 


2011; 


Madura 20ld Madura et al. 2011, 2012; Madura & Groh 2012 


; Madura et al. 


2012; Parkin et al. 


2009). 


Some other papers claim an opposite view, namelv w ~ 90° (e.s.. Abraham & Falceta-Goncalves 


2007: 


Falceta-Goncalves et al. 20051 Kashi & Soker 2008. 2009a). In Kashi & Soked (20081) we studied a varietv of 



observations of the system, and concluded that the longitude angle is w ~ 90°. Among the observations that 
led to this conclusion was the observation of the h igh-excitation narrow [Ar III] A7135 line, that originates 
in the WBs, across the 2003.5 spectroscopic event ( Damineli et al. 20081) . The high-ionization nature of the 



[Ar III] A7135 line implies that the ionizing radiation must be supplied by the hot secondary 0-star, rather 
than by the primary star. By comparing the rate of the companion's ionizing photons that reach WB D as 
a function of the orbital phase to the observed intensity of this line, we showed that a match can only be 
obtained for a longitude angle of w ~ 90°. 



Mehner et al.l (|2010l ) also presented HST/STIS spectra of high ionization lines ([Fe III] and [N III]) 
around rj Car extending over more than a full cycle. They also found that most of the emission arrives 
from the WBs. Moreover, they found a variation of the intensity of the lines over the entire orbital cycle, 
in add ition to the fa s t incr ease before periastron that followed by the typical decrease to the spectroscopic 
event. I Mehner et al.l (|2010r ) identified a fast blueshifted component of the high ionization line s that appears 



concen trated near the stars and elongated perpendicular to the bipolar axis of the system. I Mehner et al 



(|2010[ ) suggested that this component is related to the equatorial outflow and/or to dense material known 
to exist along our line of sight to the system. 

Our motivation to discuss the WBE and its interaction with the primary wind is discussed in section 
[21 In section [3] we argue that the behavior of some recently observed forbidden lines matches that expected 
from the primary wind interaction with the WBE. Our discussion and summary are in section |4l 



MOTIVATION 



Two directions motivate us to discuss the WBE and its interaction with the primary wind. First, as the 
Weigelt blob environment (WBE) moves very slowly relative to the primary, the primary wind is strongly 
shocked when it hits the WBE. The WB s are located in the eq uatorial plane in the side closer to the observer. 
In our preferred orientation of w ~ 90° (jKashi fc Sokeiil2008r ) the primary wind is constantly shocked when 
colliding with the WEB. In the oj ~ 270° case, where the secondary is toward us during most of the time, 
the secondary tenuous and faster wind hits the Weigelt blobs during most of the time. However, even in 
that case so me slower and denser primary wind gas flows toward us. More than that, the WBE extends 
tangentially ( Chesneau et al. 20051) . over an angle larger than that covered by the secondary wind after it 
collides with the primary wind. Hence the undisturbed primary wind directly collides with part of the WBE 
even in models with uj ~ 270°. We are motivated to study what are the signatures of the wind- WBE collision 
and why it is not commonly considered when modelling emission from that region. 

The second direction that motivates us comes from our view that a simple model that ignores the 
wind- WBE collision fails to account for the behavior of the [Fe HI] A4659 forbidde n line from the W BE. 
We refer to the mod el that was presented by Madura et al. ( 2011 ). and later used by Gull et al. (|201l[ ) and 
Madura et al.l (|201lh . who considered only the collision between the two stellar winds and their emission from 
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Comparing o bservations (right panel) with the mo del (left panel) propo sed by llMadura et al.ir201ll : a similar figure 
1 iMadura et"al]|2012h . This figure is taken from figure 4 of iMadura et al.l 1 I2OIII') . 



with our marking added. Presented is the 
velocity profile along a slit more or less along the equatorial plane at phase 0.976. The vertical axis is the position along the 
slit relative to the center, where up is north-east. The horizontal axis is the velocity. The different features are discussed in the 
text. 



an extended region up to ^ 1500 AU, but not the wind-WBE collision. The observations show forbidden 
line emission from a projected regio n of > 1000 AU that encloses the WBs t hat reside at real distances 
of up to ~ 850 AU from t he center (Borland et al.l 12004 : Chesneau et al.l 120051 ). assuming the WBs are in 



the equatorial plane. Both IMadura et al.l (|2012| ) and 



Gull et al 



(|201l[ ) compare their results with 3D SPH 



simulations of the stellar colliding winds (see also lGull et al.ll2009[ ). and attribute the low velocity component 
of the [Fe III] A4659 line (and some other forbidden lines) to the WBs, while the high velocity components are 
considered to result from the outflowing colliding stellar winds. They then argue that a good fit is obtained 
only for cj ~ 270°. 

We find their model to fit the observations unsatisfactory. We first examine the ob servation in wh i ch the 
spectroscopic slit is positioned more or less along the equatorial plane as presented by IMadura et al.l (j201ll ) 
at phase 0.976 (48 days before the 2003.5 periastron passage). The observations are on the right panel and 
their best model is on the left panel of Fig. [TJ In each panel the vertical axis is the position along the slit 
relative to the center, where up is north-cast. The horizontal axis is the velocity. The distance marked on 
the figure is the shock distance in our proposed model. We mark it here to emphasize that in our model the 
distance of the extended region from the center is determined by where the primary wind is shocked against 
old ejected gas. This is discussed in section [3] 

We note these discrep ancies and puzzling properties. The letters correspond to those we mark on Fig. 
|TJ (a) IMadura et al.l (|2011r) attribute the zero velocity component to the Weigelt blobs, the elliptical vertical 
white-dashed line marked by them in the right panel, while the fast components in their model come from the 
primary wind. It is not clear why two different components in their model, the primary wind and the Weigelt 
blobs, show very similar intensity, (b) It is not clear why the two different components are at about the 
same distance from the center and h ave about the same si ze, (c) There are observed redshifted components 
that their model does not reproduce. IMadura et al.l (|2012l ) discuss this in length (their appendix A; see also 
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They examine the case that the red components are contaminated. There is definitely a large 
contamination of the [Fe III] A4659 red part by the [Fe II] A4665 line. However, we note that part of the 
re d [Fe III] A4659 em ission is also seen in the lines of [Fe III] A4702 and [Nil] A57 56 in the imag es presented 
bv lGuU et all |201l[ ). and in [Fe III] A4702 also in iMadura et all |2012l ) (see also lMaduralboiOl ). Therefore 
this redshifted component seems to be real, even if much weaker than the blue shifted component and being 
heavily contaminated, (d) The observed shape of the extended emission is not as their model predicts. 
The observed extended emission is at about the same distance from the center, rather than having a larger 
distance toward lower velocities as their model predicts, (e) The observed intensity does not decrease toward 
lower velocity as their model predicts, (f) T he observed ve l ocity of the central part has a more extended 
redshifted component than a blueshifted one. iMadura et al.l (|201l[ ) model has an opposite behavior. 



The blueshifted component of the [Fe III] A4659 line is extended in the NE-SW direction and through 
the center. The zero-velocity component is extended in the same direction but to the NW of the center 



Mehner et al. 


2010; 


Gull et al. 


2011 



other in the same direction and overlapping. The redshifted component is on the NW side of the center, 
with its length as the width of the zero- velocity component. The zero velocity and blueshifted components 
seem to know about each other, i.e. they are correlated. 



Like IMadura et al.l (|201lh , iGull et al.l (j201ll) and IMadura et al.l (|2012| ) attribute the zero- velocity com- 
ponent to the WBs and the blue and red shifted components to the colliding stellar winds. As the zero 
velocity and blue shifted components seem to share some structur al features, they should have th e same 



and 


Madura et al. (2012h. I 


2011 


, 2OI2I Gull et al. 


2011) 



( 2012h. It therefore seems that the interpretation of the [Fe III] A4659 line (jMadura et al 



1 et al. (2011) 



Motivated by the considerations above we propose that the blue and redshifted components of the 
[Fe III] A4659 line (and some other forbidden lines) are formed in the postshock gas of the primary wind 
when it collides with the WBE. This explains the similar morphologies, sizes, intensities and other properties 
as we now turn to explain. The schematic proposed flow structure is drawn in the lower panel of Fig. |21 



THE PRIMARY WIND COLLISION WITH THE WBE 



3.1. The freely expanding wind 

At electron densities of < 10^ cm~'^ the emissivity ( erg s~^ cm^^^) of the [Fe III] A4659 forbidden line 
is proportional to the Fe III ions numb er density, n jjj, time s electron density, while at higher density the 
emissivity becomes proportional to n/// ( Garstang et al.iri978h . There are several quantities that are difficult 
to calculate and hence introduce uncertainties in estimating the power of the line. These include the fraction 
of the iron atoms that arc in the Fe III ionization state and the temperature of the gas. These quantities 
depend on the history of th e gas, e.g., heating in a shock and subsequent radiative cooling, and the ionizing 
radiation. iGuU et al.l (|201ll ) attempted to calculate the [Fe III] A4659 line brightness but derived values that 
are one and a half magnitude above observed values. We here limit ourselves to (a) show that our model 
can explain the power in this line, and (b) understand the morphology. We now turn to consider the central 
pi xel of size . 1" x .1" = 230 AU x 230 AU (for a distance of 2.3 kpc to 77 Car) in the observations presented 
bv lGull et al] (|201ll) . 



For an efficient emission of the [Fe III] A4659 forbidden line the iron should be in large part in the 
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The zero-velocity component 

Observation by Gull et al. (2011) 




View from the direction marked by the 
arrow (to the positive y direction) 



Observer 



The blue-shift component 





The WBE 



Shock wave 



Primary 
wind 



The red-shift component 



The binary system 



Fig. 2. — Upper panel: The [Fe III] A4659 forbidden line from lGuU et al.l ll201ll'l . The emission is presented in three velocity 
components: blueshifted component (—400 to —200 km s~^), zero-velocity component (—90 to —30 km s~^) and red-shifted 
component (-1-100 to +200 km s~^). Bottom panel: View from the direction marked by the white arrow, according to our 
model. In this model the emission comes from the shocked primary wind, that collides with the WBE — the dense slowly 
moving gas near the equatorial material. The flow direction is schematically drawn, as well as the location of the observer 
according to our model. The axes are defined in the figure: The x-y plane is the projected view on the sky where the x-axis is 
more or less along the projected axis of the Homunculus on the sky. The z-axis is directed to the observer. 



ionization state of Fe I II, which it turn requires a temperature of several xlO'^ K ^ T < 10^ K (e.g, 
Gnat &: Sternbersll2007l ) . The wind is not at an exact equihbrium as it cools during expansion. However, 
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the regions of the wind that are relevant to us are those that are ionized and heated by the secondary 
star. We assume that in these regions the coohng is not much faster that the recombination and ionization 
time scales of the relevant ions, and hence can be considered for our purposes as being in equilibrium. In 
our model the basic process is that the emitting gas is at ~ 10"* ^ K and more or less in equilibrium. At 
that temperature most of the iron is in the Fe II ionization state (jGnat fc Sternberg] I20Q7I ). The secondary 
radiation then ionizes the Fe"'""'^ ions to form the desired Fe"*"^ ions. 

The central pixel has a size of 230 AU, reaching an average distance of 130 AU from the center, which is 
much larger than the wind interaction zone. We therefore consider the density of the freely streaming primary 
wind (for the colliding zone contribution see below) . The LBV primary star blows a wind with a high mass 
loss rate of Mi — 3 x 10~ '* Mq yi"-^ and with a term i nal ve locit y at the equato r of vi c=l 500 km s~^ (e.g., 
Pittard &: Corcoran 2002). We note that Gull et al. ( 2011 ) and Madura et al. ( 20121) take a higher mass 
loss rate of ~ 10~^ Mq yr~^. However, their SPH simulations lack radiative cooling, and therefore their 
postshocked densities are underestimated. We consider distances much larger than the orbital separation and 
neglect the motion of the primary around the center of mass. The electron density of the freely streaming 
primary wind is 



riei 7 X 10^ 



Ml 



3 X 10-4 M© yr- 



(lOOAu) ( 



Vl 



500 km s- 



(1) 



We now scale quantities to estimate the power of the central pixel. In planetary nebulae the power 
of the [Fe III] A4659 line reac hes values of < 0.025 times the power (e.g., Delgado Inglada et al. 2009t 
Simon-Diaz fc Stasinskal 120111 ). We checked and found that such a ratio will be too low to explain the 
observed intensity of the [Fe III] A4659 line in ij Car. This is explained by the small rat io of Fe III ion 
numb er-density to hydrogen number density in planetary nebulae nm /uh ^ 3 x 10~ ^ ( Garstang et al 



19781 ). In other environments the fraction can be much higher, e.g., ICarstang et al.l (|l978l ) find this ratio to 
be ~ 3 X 10"^ in the Orion nebula. We scale with this ratio, which implies that ~ 10% of the iron is in the 
Fe III ionization state. To reach this fraction the hard ionizing radiation of the seconda ry sta r is required. 



The emissivity in the line at densities of rip. ^ 1 0** cm ^ according to Garstang et al. (|l978t we note the 



updated calculations of, e.g. jKeenan et al.lll992l but for the purpose of our estimates the simple expressions 
of Garstang et al. are adequate) is 



em ^ 3 X lQ-''piii{T) 



niii 



3 X lG-6nif / VlO^ cm 



erg s cm 



(2) 



where 6 — 2 ior lower densities and S ^ 1 for higher densities. The coefficient f3jjj{T) has values of 3, 1, 0.06 
for temperatures of 2, 1,0.5 x lO'' K, respectively. 

Equations ([Ij and (0) imply that outside the central pixel, where (5 ~ 2, the surface brightness of the line 
decreases with radius according to 5 oc r^^, and the central pixel will be much brighter than its immediate 
surroundings. We perform the integration over the emissivity from radius Tmin to the average distance of 
the first pixel 130 AU, and take 6 = 1.5 as a gross average in that range. Taking a distance of 2.3 kpc to f] 
Car we find the flux from the central pixel to be 



10-"/3///(r) 



Nhi 



Ml 



1.5 



3 X IQ-f^riH 



500 km s- 



-1.5 



3 X 10-4 Mo yr-i 
130 AU 



(3) 



erg s ^ cm ^. 
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We no te the dependence on the primary wind parameters. For a mass loss as high as used by I Gull et al 
(j2011[ ). for example, the flux will be 6 times higher than the estimate here. 



Had we considered the coUiding-wind region emission would have been larger, both due to the higher 
density of the shocked primary wind and due to the strong ionizing radiation that can increase the density 
of Fe"*"^ ions. The reason fir this is that the higher density implies also that the dependence on density 
(equation [5]) is linear {6 — 1), rather than having the steeper dependence oi 6 = 2. Though the density of 
the shocked primary wind is higher, the volume occupied by it is smaller than if it was freely streaming. 
Therefore, we expect the addition from the shocked primary wind to be only slightly more had it been freely 
streaming. 



The observed flux in the central pixel observed bv lGull et al.l (j2011l) is ~ 3 x 10^^^ erg s^^ cm^^ (when 
summing up the three velocity components). Despite the crude derivat ion we can conclude that the primary 
wind can account for the brightness of the central pixel as observed bv lGuU et al.l (|201lh . The primary wind 
explanation accounts also for the central pixel being much bri ghter than its i mmed iate environment (not 
including the extended bright regions). The model presented bv lMadura et al.l (|2012r) and I Gull et al.l (|201lh 
does not reproduce a bright central region that is detached from the extended emission. 



3.2. The Weigelt Blob Environment (WBE) 



Consid er WB D as an exa mple. By inspecting the relevant figures in lDorland et al.l (|2004l ). IChesneau et al 

( 2005 ) and Gull et al. ( 20091 ). it can be seen that WB D covers an opening angle of aw — 30° of the 
primary wind. The radius of WB D is rf, ^ 220 AU, and its distance from the center is taken to be 
~ 850 AU. The electron de nsity in WB D is c omparable to the density in the other blobs and estimated to 
be nw ^ 10^ - 10^" cm'^ jVerner et al.ll2002f ). The mass of WB D is Mw ^7 x 10-^(71^^/10'^ cm'^) Mq. 
Later we note that for our model the much extended, somewhat lower-density region around WB D plays 
the main role in colliding with the primary wind. 



From HST/STIS spectra of Fe II, [Fe II], [Ni II] and [Ni III] lines from WBs C and D. ISmith et al.l (|2004 
measured the radial velocities of these blobs, and found them to be vw ^ 40 km s^-^. With an opening angle 
of aw = 30° the blob intersects a fraction of fw = 0.017 of the primary wind. Over the 120 years since the 
Lesser Eruption this a mounts to a mass o f AMw — 6 x 10^^(/u)/0.017) Mq from the primary wind that 
collided with the blob. Raea et al. ( 19981 ) have shown that interaction of a wind with a static molecular 
clump can cause the clump material to accelerate without being dissociated, thus producing high-velocity 
molecular emission. It is therefore safe to assume that the interaction of the primary wind with the blobs 
can accelerate them. The primary wind would accelerate the blob to a speed of Va — vi{AMw /Mw)- The 
mass of the blob should therefore be Mw ^ 0.01 M fr^ in order for it not to be accelerated to fa > 40 km s^^. 
The WBs are part of a more massive gas complex (jChesneau et al.ll2005l ). which we term the Weigelt blob 
environment (WBE), that withstands the acceleration by the primary wind. 

As the WBE was not accelerated to a high velocity, we can put a lower constraint on its mean density 



J^WBE,!!!!!! 3 X 10 



10 



Ml 



3 X 10-4 Mq yr- 



- 200 AU 



(4) 



We estimate the density of the central part of the WBE to be nwBE ~ 10^ cm-'^. In the outskirts of the 
WBE the density will be lower than nwBE , due to interaction with the primary wind and evapor a tion. The 
inner edge of the WBE is closer to the binary system. From images given by IChesneau et al.l (j2005l ) we 
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estimate the average value of the distance between the inner edge of the WBE and the binary system to be 
about half the distance of the center of the WBs or less, rd < 400 AU. 



3.3. Interaction with the WBE 

To check whether the interaction between the primary wind and the WBE can produce the observed 
flux from the [Fe III] A4659 line, a simple estimate can be done as follows. The inner edge of the WBE is at 
a distance of Vd from the binary system. The immediate post shock electron density of the primary wind is 
^ 4 times the pre-shock electron density, and it is given by 

^ 2 X 10^ f -\ f ' { ^) (5) 

^ \i X 10-4 M© yr~i j V400 AU/ V500 km s-^ ) ^ ' 

The shocked gas is heated to a temperature of up to T ~ 3 x 10^. It cools radiatively to 10^ K within a 
time scale of ~ 6 months. However, when it starts to cool it is compressed, and the cooling time shortens. 
Overall, the gas cools over a time scale of a few months. The flow time is i/ > r^/wi ~ 4 yr, such that the 
post-shock gas cools quite rapidly. The gas cools to ^ 10'* K as it flows. This is the time when it efficiently 
emits the [Fe III] A4659 forbidden line if it is ionized by the UV radiation from the companion star. A 
cooling by a factor of ^ 100 can result in a compression by a similar factor near the stagnation point of 
the interaction. However, the gas flows outward from the interaction region and the density decreases. We 
scale the density with an increase by about one order of magnitude from the immediate post shock region 
to Uef ^ 10^ cm~^. 

The size of an arcsec^ at a distance of 2.3 kpc to -q Car is (2300 AU)^. We take the length of the region 
to be ~ Td- By using equation ([2]) we find the surface brightness of the postshock primary wind in the 
[Fe HI] A4659 line to be 



The maximum flux in the [ Fe HI] A4659 l i ne fro m a pixel of 0.1" x 0.1" as presented bv iGull et al 

(|201l[ ) is ^ 10^*^ erg cm"'^ s^^. iMehner et al.l (|201ll ) find the maximum power to be from WB C, with a 



flux of 2 X 10~^^ erg cm^^ s^^ within 0.1" from the blob. Our simple estimate shows that the postshock 
primary wind can account for the intensity of the [Fe HI] A4659 forbidden line in the extended region within 
one arcsec from the center. We emphasize that the postshock region in our model for the extended emission 
refers to the shock formed when the primary wind collides with the WBE and not with the secondary wind. 



3.4. Geometry 

We consider the 3D geometry of the primary wind interaction with the WBE, and take i nto account the 



inclination angle and the orbital orientation. The WBs extend towards us and to the NW side (jPorland et al 



2004 IChesneau erZlboOSl: ICufl et al.llioogi) . In Fig. [5] and [3] we present the geometry of the WBE. We 



use an axes system where the x-y plane is the plane of the sky, where the x-axis is more or less along the 
projected axis of the Homunculus on the s ky, and the z-axis is directed to the observer. The observer in 



our model, as in our previous papers (e.g., iKashi fc Sokeiil2008l l2009at l is located at a longitude angle of 



u) — 90°, namely, the primary is towards the observer for most of the orbit. The secondary is towards the 
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observer only for a short time cfose to periastron passage. As the WBE is located between the observer 
and the binary system, the primary is also towards the WBE for most of the orbit. Therefore, the WBE is 
exposed to the primary wind which collides with it, goes through a shock wave, and flows around the WBE. 



In our model the emission of the [Fe III] A4659 forbidden line, as well as some other lines such as the 
[Fe III] A4702 and the [Nil] A5756 lines, come from the shocked primary wind, after it collides with the WBE. 
One important note is in place here: the reader should be aware that usually the term 'shocked primary 
wind' refers to the primary wind after it is shocked when colliding with the secondary wind. But here we 
refer to the shocked primary wind after it collide s with the WBE. The flow direction of the wind after this 
collision is schematically drawn in Fig. [5] and |31 I Gull et al.l (|2011l ) divide the [Fe III] A4659 forbidden line 
emission into three components: blueshifted component (—400 to —200 k m s~^), zero - veloc ity component 



(-90 to +30 km s"^) and red-shifted component (-1-100 to -1-200 km s"^). iGull et al.l (|2011l ) attribute the 
zero velocity component to WBs. We accept this in our model, and more specifically we attribute the zero 
velocity component to the edge of the WBE region which is impinged by the primary wind and ionized by 
the secondary ionizing radiation. The WBs themselves are too dense to emit the [Fe III] A4659 forbidden 
line, but the density at the WBE edges is low enough to allow it. We can see this from the constraint derived 



above (Section 13. 2p on the density of the WBs which is 
efficient. 



10 cm , where the emission of the line is not 



The blueshifted and redshifted components come from the postshocked primary wind gas after it cools 
to ~ 10** K and is ionized by the secondary star. The division into the blue and red shifted components is 
caused by the postshock flow on the two sides of the equatorial plane of the binary system, where the WBE 
resides. Because the line of sight is inclined to the binary axis the two components are not symmetric in the 
position- velocity diagram. As shown in equation ^ and the related discussion, this region can quantitatively 
explain the flux in the emission line. 



The peak radial velocity component of the forbidden emission lines is at ^ 400 km s ^ (jCull et al 



201 II ). lower than the terminal velocity of the primary wind of ^ 500 km s~^. In our model the slowing down 
is a consequence of the interaction of the primary wind with the slowly moving WBE. In a model where 
the winds are attributed to the primary wind collision with the faster secondary wind, higher velocity than 
the primary wind velocity might be obtained. Indeed, the model of Gull et al. ( 201l[ ) predict blueshifted 
velocities in the extended components that are ~ 20% higher than observed (see flgure[T]). 



DISCUSSION AND SUMMARY 



As the primary stellar wind collides with the slow dense material in the Weigelt blob environment 
(WBE) it is shocked. The postshocked gas reaches a temperature of ~ 10^ K but rapidly cools to ~ 10** K 
and is compressed. This warm dense gas is expected to become a source of some emission lines, including 
forbidden lines. The interaction of the primary wind with the WBE cannot be ig nored whe n mod elling 
emission from that region. The high resolution observations reported recently by iGull et al.l (|201l[ ) and 



Madura et al 



(2012) that extend to ~ 0."5 from the binary system overlap with the WBE as mapped by 



Chesneau et al.l ( 20051 ). We showed that th e intensity and loca t fon of the differe nt velocity components of 



the [Fe III] A4659 extended forbidden line (jMadura et al.l 1201 ll : iGuU et al.ll2011r) . can be explained by the 
interaction of the primary wind with the WBE. The geometry of our proposed model is depicted in figures 
[5] and 121 while the intensity is calculated in section [ 
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Observer 




L _^>qutarSiM_( 1 ^ Va(ea(y(lanB'^) 

Fig. 3. — The geometry of the WBE and the source of the [Fe III] A4659 forbidden line emission. The bottom panel, taken 
from iMadura et al shows observations of the line in 3 position angles. WBs B, C, and D are marked as black points. 

The upper panel and the arrows connecting the bottom panel to the upper one show from where each of the three components 
is emitted according to our model. The zero-velocity component comes from the edge of the WBE region which is impinged 
by the primary wind and ionized by the secondary ionizing radiation. The blueshifted and redshifted components come from 
the postshocke d primary wind gas after it cools to ~ 10^ K and ionized by the secondary star. The redshifted component is 
contaminated jMadura et al.ir2012i) . However, as discussed in the text, we find that a weak redshifted component does exist on 
top of the contamination. 

In our proposed model the freely expanding primary wind continuously collides with the WBE and is 
shocked. As the WBs are closer to the observer, so should be the primary most of the time. This means 
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that the orientation is such that lo ~ 90°: the primary is closer to us for most of the binary orbit. The 
secondary becomes closer to us only for a short time near peria s tron p assage. This model overcomes some 
difficulties in the model of iMadura et al.l (|2011l ) and iGuU et al.l (j201lh that assumes an opposite longitude 
angle of w ~ 270° where it is the secondary that faces the WBE for most of the time (see section [5]). 



Our explanation of the recent observations of iMadura et al.l (120111 ) and lGull et al.l (|201ll ) with a model 
where the primary star is closer to us for most of the orbit adds to a growing number of observations that 
support this orientation. The observations that support an orientation of a; ~ 90° can be summarized as 
follows. 

(1) The Doppler periodicity of the P Cygni component of the He I optical lines is well explained if the origin 
of the absorption is taken to be in the acceleration zone of the secondary wind (jKashi fc Sokeiil2008l ). and 
the binary system is oriented such that lu ~ 90°. I.e., the secondary wind absorbs from the emission of the 
secondary photosphere. 



(2) The Doppler variation of the N II AA5668-5712 lines observed bv lMehner et al. ( 2011 ) i s also well ex- 



plaine d if they come from the acceleration zone of the secondary wind ( Kashi fc Sokei] 2011 ) . Mehner et al 



(|201lt ) criticized this claim by the fact that the very small locale near the secondary star cannot produce 
substantial absorption lines i n the spectrum of the much more luminous primary. Our answer to that criti- 
cism was already discussed in Kashi fc Soker ( 20081 ) in the context of the He I lines. The argument goes as 
follows. The absorption in both the N II and the He I lines is < 15% , and in some cases is from the excess 
emission of the P Cygni lines (see observations in lHillier et al.ll200ll ). The secondary contributes ~ 20% of 
the luminosity in the optical wavelengths and therefore there is no problem for the acceleration zone of the 
secondary to account for the absorption intensity of these lines. 

(3) The Doppler shift of the low ionization Fe II A6455 line is well explained by a a; ~ 90° orientation when 
it is attributed to the primary wind (Kashi fc Soker 2o"oi). 

(4) The periodic Doppler shift variation of the hydrogen Paschen lines is also explained by the same orien- 
tatio n if the lines are attributed to the apex region (stagnation point) of the colliding winds (jKashi fc Soker 
20081) . 

(5) The observed variation of the He I A10830 line ( Groh et al. 2010l) . in which the f ast blue absorption 
component exists for only several weeks prior to the periastron passage was explained in lKashi et al.l (j2011l ) 
using a 3D numerical simulation. This simulation showed that the fast variation is well reproduced for 
w ~ 90°. This model took into account the fact that the line is emitted from an extended source around 



the primary, an ingredient that was missing in the model of iGroh et al.l (j2010[ ) who deduced an opposite 
orientation. 

(6) T he hydrogen column density deduced from X-ray absorptio n in XMM - Newton observations (iHama guchi et al 
20071 ) could only be explained quantitatively by this orientation (jKashi fc Soker 2009a l . It was quantitatively 
shown that an orientations of cj = 0°, 180° and 270° give results contradicting observations. It is important 
to note that in contrast to the hydrogen column density, the X-ray light curve by itself cannot be used to 
deduce the periastron longitude. Minor adjustments of some parameters of the binary system (inclination, 
eccentricity, etc) allow one to fit the X-ray light curve with almost any orientation (see further discussion in 
Kashi fc Sokei]l2009al) . 

(8) T he variation in the in tensity of the highly excited [Ar HI] narrow lines, which also originated in the 
WBs ( Damineli et al. 20081) . follows the ionizing radiation from the companion for an orient ation where the 
second ary is towards the WBs only for a short time at periastron passage. In section 4 of iKashi fc Soker 
(j2008l ) we showed that even when behind the primary the secondary can ionize the WBE. The reason is 
that the primary w ind absorbs the high-ionizing secondary radiation only within a small angle close to the 
primary (jKashi fc Soker, ,2007i) . and the WBE is extended. The fact that in our model the primary wind 
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absorbs a non- negligible (but not all) o f the secondary ioniz ing radiation, allows us to explain the varia- 
tion within a cyc le, a nd between cyc les ( Kashi fc Soken 120081 ) . On the contrary, in the model preferred by 
Gull et al.l (|2011r) and iMadura et al.l (|2012r) . the secondary faces the Weigelt blobs. This implies that the 
Weigelt blobs received a constant ionizing emission for most of the orbital period, hence a constant [Ar III] 
emission is expected. This contradicts observations. 

The determination of the orbital orientation, i.e., whether the primary or the secondary is closer to 
us at periastron passage, has implications that go beyond the goal of explaining the periodic behavior 
of ry Car. Knowledge on the orientation will allow us to better understand the interaction of the binary 
system near periastron passages. The emission and absorption behavior of 77 Car near periastron passage 
cannot be explained if the secondary wind is not substantially suppressed then. This suppression of the 
secondary wind is thought t o be caused by accre t ion of primary wi nd material onto the secondary star 
near periastron passage fe.g.. lKashi fc Sokeiil2009bt lAkashi et al.ll201lh . Numerical simulations show indeed 
that near periastron passage instabilities in t he wind colliding r egion le ad to the formation of dense blobs 
that are accreted toward the secondary star (jAkashi et aLluOllt but see lParkin et al.ll2011l ). Although the 
requirement for the accretion process does not depend on the orientations, knowing the orientation will help 
to understand the onset of, and the exit from, the several week-long accretion phase. Understanding the 
accretion process in present-day i] Car might shed light on the major accretion process that took place during 
the 1837-1856 Great Eruption, with implications to other LBV binary systems. 

We thank one of the two referees (A. Moffat) for very helpful comments that substantially improved the 
manuscript. This research was supported by the Asher Fund for Space Research at the Technion, and the 
Israel Science Foundation. 
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